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ABSTRACT  . 

Each  of  the  two  principal  methodologies  for  energy  analysis 
has  drawbacks  which  may  severely  limit  the  accuracy  of  results.  A 
method  for  combining  the  two  methods  is  presented,  and  is  shown  to 
yield  results  more  accurate  than  either  method  used  independently. 


INTRODUCTION 

This  paper  discusses  two  methods  for  calculating  the  energy- 
cost  of  individual  goods  and  services,  or  more  complex  systems.   It 
presents  a  method  for  combining  them  in  a  way  that  minimizes  uncer- 
tainty in  the  result. 

The  "energy  cost"  of  an  item  is  defined  as  the  total  energy 
required,  directly  and  indirectly,  to  produce  it.   For  example,  the 
energy  cost  of  a  beer  can  includes  not  only  the  energy  consumed  by 
the  can  manufacturer,  but  also  that  required  to  smelt  the  aluminum, 
to  mine  the  bauxite,  plus  all  the  energy  needed  to  transport  these 
intermediate  products. 

Energy  cost  is  usually  calculated  using  one  of  two  methods. 
The  first,  called  process  analysis,  starts  with  the  final  product 
and  identifies  direct  process  inputs,  then  the  inputs  to  those  products, 
etc.   This  "tree"  of  process  inputs  is  shown  in  Fig.  1.   The  energy 
cost  is  obtained  as  the  sum  of  the  direct  energy  inputs  at  each 
juncture.* 

The  other  method,  based  on  input-output  analysis  starts  with 
a  description  of  an  N-sector  economic  system  producing  N  types  of  goods 
and  services.   The  energy  cost  of  all  are  determined  simultaneously  as 
the  solution  of  a  set  of  N  linear  equations.** 


*The  direct  energy  input  to  a  process  is  the  enthalpy  of  the  fuel  con- 
sumed in  that  process.  For  examples  of  process  analysis  see  Hannon  (1973) 

**The  adaptation  of  economic  input-output  theory  to  energy  analyses  is 
described  by  Bullard  and  Herendeen  (1975). 
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Figure  1.   The  PROCESS  ANALYSIS  "TREE 
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The  data  required  for  each  method  are  identical;  the  tech- 
nology* for  producing  all  types  of  goods  and  services  must  he 
specified.   The  energy  cost  per  unit  of  output  in  each  sector  or 
energy- intensity  is  a  function  of  this  specification  alone. 

Process  analyses  are  tedious  and  are  usually  truncated  after 
only  a  few  steps.   This  is  usually  done  to  limit  data  requirements  to 
those  describing  production  of  only  the  most  important  inputs  to  a 
process.   Truncation  errors  are  unknown.   Input-output  (1-0)  analysis 
avoids  the  truncation  problem,  as  it  includes  all  levels  of  inputs 
and  feedback  loops.   It  is  limited,  however,  by  the  degree  of  dis- 
aggregation of  the  model  which  is  necessarily  general  and  not  problem- 
specific.   For  example,  1-0  gives  the  energy  cost  of  metal  cans  -  not 
beverage  cans;  of  motor  vehicles  -  not  Chevrolets . 

The  purpose  of  this  paper  is  to  present  a  procedure  for 
combining  the  two  methods  to  minimize  the  disadvantages  of  each.   The 
combination  is  straightforward  because  both  methods  are  linear  and 
both  require  the  same  type  of. input  data.   Finally,  this  "hybrid" 
analysis  technique,  to  systematically  reduce  uncertainty  in  the  energy 
cost  of  a  commodity,  is  applied  for  an  example  calculation. 

The  "Ideal"  Process  Analysis 

Consider  the  most  general  case  of  an  economic  system  in  which 
there  are  N  types  of  goods  and  services  (N  may  be  arbitrarily  large, 


*In  this  paper  the  word  technology  is  to  be  interpreted  in  a  narrow  sense, 
as  describing  material  and  energy  inputs  to  production  processes.  A  more 
detailed  definition  is  given  in  the  next  section. 


-  k  - 

on  the  order  of  thousands  or  millions).   In  order  to  calculate  the 
energy  cost  of  one  of  these,  a  diagram  or  tree  is  constructed  such 
as  that  in  Fig.  1.   There  may  be  up  to  N  inputs  at  each  node  and  the 
number  of  nodes  is,  in  principle,  unlimited.  At  each  node  the  direct 
energy  input  is  tabulated;  the  energy  cost  of  the  good  or  service  is 
the  sum  of  these  inputs. 

For  a  typical  product,  n,  the  production  technology  is  repre- 
sented by  a  vector  a  where  a  typical  element  a.   represents  the 
amount  of  product  i  needed  directly  to  produce  a  unit  of  product  n. 
The  N  x  N  matrix  a  then  provides  a  linear  representation  of  the 
technology  of  producing  all  goods  and  services. 


by: 


Let  e  represent  energy  intensity  of  product  n.   It  is  given 


N  N   N 

e  =6   +a   +  I  a.  a.  +  Z   E  a  „  a. .  a.  +••••    (l) 

n    en    en   ._   ei  in   ,  ,  .  ,   ej  ji  in 


where  the  subscript  e  denotes  the  energy  sector  and  6   =0  for  n#e; 

6   represents  the  heat  content  of  a  unit  of  energy.*  The  term  a 

ee  *  °*  en 

denotes  the  energy  used  directly  in  producing  a  unit  of  product  n,  and 
succeeding  terms  correspond  to  direct  energy  inputs  at  each  level  of 
the  process  analysis  tree. 


*Energy  is  usually  measured  in  terms  of  enthalpy,  so  6   =  1. 
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In  practice,  such  a  large  number  of  terms  is  never  computed. 
At  the  first  level  only  the  most  significant  inputs  are  considered, 
and,  of  those,  only  a  subset  is  further  broken  down  into  its  components, 
Unfortunately,  diminishing  contributions  from  each  level  provide  no 
information  for  a  truncation  rule.* 

Input-Output  Analysis 

This  is  a  technique  for  representing  the  entire  system  of 
N  production  processes  as  a  linear  network  model.   Each  node  charac- 
terizes a  sector  of  the  economic  system,  each  producing  a  unique  good 
or  service.  Figure  2  shows  the  energy  flows  entering  and  leaving  each 
sector. 


N 

E  e.  T. 

.  ,   1  in 
i=l 


I 


E 


n 


Figure  2.   Energy  Balance  for  a 
Producing  Sector 


Energy  embodied  in  inputs  from  other  sectors,  e.T.  ,  enters 
at  the  left,  while  energy  embodied  in  the  sector's  output  X  is  shown 


^Consider  for  example,  the  divergent  series  E— . 
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exiting  at  the  right.   If,  in  Fig.  2,  sector  n  denotes  the  energy 
sector,  an  amount  E  is  extracted  from  the  earth.   The  energy  "balance 
equation  becomes: 


N 

I     e.  T.   +  E  =  e  X  (2) 

,  ,  .i  in   n   n  n 

i=l 


or,    in  matrix  notation  we  have: 

Ll    +     I    =     e.i  ,  (3) 

which  is  a  set  of  N  equations  that  can  be  solved  for  the  N  unknowns  e_. 
2£  is  the  diagonal  matrix  whose  elements  represent  the  total  output  from 
each  sector. 

From  the  definition  of  ^  given  earlier,  we  have: 

and  eq.  (3)  becomes: 

L=    £.  (I-^)"1  (5) 

where  e  is  a  unit  vector  which  identifies  the  energy  sector  row  of 
(L-ji)   as  the  energy  intensities. 
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Combining  Process   and  Input-Output  Analyses 

It  is  now  clear  that  the  energy  cost  of  any  good  or  service 
can  be  determined  by  either  process   analysis    (eq.    l)    or  input-output 
analysis    (eq.    5).     Both  methods  are  linear  and  require  identical  input 
data.      The  proof  that  the  results  are  identical  is   straightforward. 

If  the  spectral  radius  of  £  is  less  than  unity,*  we  have 

£   =   £     [£  +  4  +  &    +  &3  + ]  (6) 

which  shows  that  the  energy  intensities  calculated  using  eq.    (5)   are 
identical  to  those  obtained  through  the  more  tedious  process  analysis 
eq.    (1). 

For  most  applications,   however,  this   rather  complete  set   of 
data   (the  N  x  N  matrix  &)   is  not   available  at  the  necessary  level  of 
detail.      It  exists  only  at  a  more  aggregated  level  K<N,  where  K~360  for 
the  United  States   economic  system,   and  is  much  smaller  for  most   other 
nations. 

Because  of  this  lack  of  data,   input-output  results  give  only 
the  average  energy  intensity  of  a  sector's   output.      Accuracy  is  limited 
by  the  level  of  aggregation:     the  energy  intensity  of  "metal  cans" 
would  apply  to  both  aluminum  beer  cans   and  55-gallon  oil  drums. 

Process   analysis   does  provide   a  framework  for  obtaining  new 
detailed  problem-specific  data.      For  example,   the  energy   consumed 


*This   condition  is  indeed  satisfied  for  the  360-sector  description 
of  the  U.S.    economic  system. 
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directly  by  the  beer  can  manufacturer  is  the  second  term  in  eq.    (l). 
However,    the  process   analysis  must  soon  be  terminated  as  the  data 
acquisition  effort  rapidly  increases   and  certain  data  are  found  to 
be  unavailable. 

Fortunately,  the  truncation  error  can  be  minimized  using 
the  results   of  input-output   analysis.      Note  that  the  series   expansion 
in  eq.    (l)    can  be  truncated  at   any  level. 

e.=  e[l  +  a=+a2  +   •••   +  a.m  (i-a)"1]    .  (7) 

It  is  trivial  to  prove  the  matrix  product  in  the  last  term  commutative, 
so  the  energy  intensities  are  given  by: 

£=e[l+a+a2+   ••••   +  a*1"1]  +  e/af  (8) 

where  e/  is  a  vector  of  energy  intensities  described  in  more  detail 
below.  To  compute  the  energy  cost  of  a  particular  item,  one  may 
evaluate  the  first  few  terms  from  available  (problem-specific)  data 
and,  in  doing  so,  truncate  the  expansion  early.  Ultimately,  the  entire 
matrix  |i  would  be  filled  as  the  number  of  steps  m  increases.   It  is 
filled  columnwise,  as  the  technology  for  producing  each  input  is 
specified.   Typically,  the  process  is  terminated  at  m=2  or  3,  before 
too  many  inputs  are  involved. 

The  last  term  in  eq.  (8)  approximates  the  sum  of  the  truncated 
terms.   The  vector  e_'  is  obtained  from  an  input-output  analysis  of  the 
economic  system  described  at  the  K-sector  level  of  detail,  where  K<N. 
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e/  is  an  N-order  vector  containing  only  K  distinct  values.   The  order 
N,  required  for  eq.  (8),  is  obtained  "by  repeating  each  of  the  K  distinct 
energy  intensities  a  number  of  times.   The  error  involved  in  this  trunca- 
tion  depends  on  the  uncertainty  introduced  by  characterizing  these  m 
level  inputs  as  "typical"  outputs  of  industries  aggregated  to  the 
K-sector  level. 

In  summary,  process  analysis  provides  a  framework  for  utilizing 
a  limited  amount  of  problem-specific  data  to  reduce  the  "aggregation 
error"  inherent  in  input-output  results.   Using  1-0  results  to  truncate 
the  analysis  eliminates  the  problem  of  an  unknown  truncation  error, 
replacing  it  by  a  smaller  aggregation  error  associated  with  energy- 
costing  the  higher  order  inputs.   We  shall  call  the  combination  of  these 
techniques  "hybrid  analysis"  and  describe  the  procedures  below. 

Error  Reduction  Criterion 

In  practice,  each  step  in  a  process  analysis  may  be  viewed  as 
an  expansion  of  the  system  boundary  (around  the  item  being  analyzed) 
into  the  economic  system,  tabulating  direct  energy  inputs  at  each  step 
(see  Fig.  3).   The  results  of  input-output  analysis  may  be  used  to 
estimate  the  energy  embodied  in  flows  crossing  the  system  boundary  at 
any  level,  by  associating  each  good  or  service  with  one  of  the  K<N 
sectors  of  the  1-0  model. 

The  1-0  results  are  indifferent  to  the  location  of  the  system 
boundary  interface.  Regardless  of  the  number  of  process  analysis  steps 
taken,  the  system  boundary  still  looks  the  same  from  the  1-0  side;   only 
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K  types  of  goods  and  services  cross  the  boundary.   The  energy  cost  of 
each  of  those  is  known. 

For  industrial  sector  i,  e".    is  the  energy  intensity  (eq.  8) 
and  a.  is  the  amount  of  product  i  crossing  the  boundary.   The  energy 
embodied  in  the  flow  of  sector  i  goods  is  therefore  e^a  .   This  flow 
term  is  analogous  to  the  truncation  term  in  eq.  (8)  for  energy  intensities 
In  general,  eT  and  a.  have  uncertainties  AeT  and  Aa.,  respectively.* 

The  criterion  for  optimizing  the  hybrid  analysis  procedure 
is  to  minimize  uncertainty  in  total  energy  costs.   For  this  purpose, 
inputs  from  the  K-sector  economic  system  at  each  process  analysis  step 
are  classified  into  one  of  three  categories:   Type  1;  typical  of  a 
sector's  output,  Type  2;  atypical  of  a  sector's  output,  or  Type  3; 
miscellaneous  (not  specified  as  output  from  a  particular  sector).  For 
atypical  inputs,  the  analyst  must  increase  energy  intensity  uncertainties 
to  AeT>AeT,  reflecting  the  aggregation  error.  Type  3  inputs  are  assigned 
an  average  energy  intensity  designated  as  e'  with  uncertainty  Ae'.** 

Obviously,  the  total  uncertainty  (e'Aa  +  Ae'a  +  Ae'Aa)  can 
be  reduced  by  reducing  the  Aa  term.   Uncertainties  in  energy  intensities 
(e')  of  Type  2  or  3  inputs  can  usually  be  reduced  by  a  process  analysis 
of  their  inputs.***  Figure  k   shows  the  first-order  errors  associated 
with  a  Type  2  or  3  input  disaggregation.   A  step-by-step  procedure  to 


*For  a  discussion  of  the  Aef  uncertainties  see  Bullard  and  Sebald  (1975). 

**Usually,  the  value  used  here  is  that  of  the  energy /GNP  ratio  if  inputs 
are  measured  in  dollars. 

***Because  Aa's  may  increase  as  disaggregation  continues,  minimum  uncer- 
tainty will  not  necessarily  be  achieved  by  greater  disaggregation. 
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reduce  the  overall  uncertainty  is  demonstrated  "by  a  simplified  example 
in  the  Appendix. 

SUMMARY 

To  reduce  the  uncertainty  that  is  unavoidable  in  either  of  the 
conventional  energy  analysis  methodologies,  the  methods  can  be  combined 
to  minimize  the  drawbacks  of  each.   Process  analysis  alone  is  intractable 
because  of  the  massive  data  requirements  needed  for  completion;  a  truncated 
process  analysis  has  an  unknown  error.   Input-Output  analysis,  on  the 
other  hand,  is  blessed  with  a  global  data  base  but  is  severely  constrained 
by  aggregation  problems  and  an  outdated  description  of  technology. 

Using  the  combined  "hybrid"  analysis  technique  presented  here, 
errors  can  be  minimized  and  quantified.   Depending  on  the  data 
obtainable,  the  complexity  of  the  analysis  can  be  traded  off  against 
the  accuracy  desired. 
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APPENDIX:   A  THREE  SECTOR  EXAMPLE 

This  section  illustrates  the  hybrid  analysis  procedure  for 
energy  costing  "widgets"  using  a  hypothetical  three-sector  economy. 
The  table  below  gives  the  necessary  1-0  data  for  the  analysis. 
Figure  5  shows  the  dollar  flows  and  system  boundary  definitions  for 
each  level  of  the  analysis.   (For  simplicity,  only  positive  errors 
are  considered. ) 


INPUT-OUTPUT  DATA  FOR  A  HYPOTHETICAL 
THREE  SECTOR  ECONOMY 


Sector 

Energy  Intensity  (Btu/$) 

Uncertainty  (Btu/$) 

1 

k.O 

.5 

2 

T.O 

1.0 

3 

15.0 

1.0 

Miscellaneous 

8.0 

10.0 

Only  positive  error  terms  are  considered  here, 
Corresponds  to  Type  3  inputs. 


Step  I.   The  initial  step  evaluates  the  energy  costs  and 
uncertainties  at  the  most  aggregated  level  of  analysis.  At  this  level 
inputs  are  all  considered  to  be  of  Type  3.  The  total  dollar  costs  and 
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un certainties  are  found  (e.g.,  from  a  widget  expert)  to  be  as  follows: 

a   =  10.0  Aa   =  1.0 

o  o 

The  resulting   energy  costs   and  associated  first-order   (positive)    error 
becomes: 

Energy  Cost   =  Z   e  a  =    (8.0)(l0.0)   =  80.0  Btu 

Positive  Error  Term  =     Z(eAa  +  Aea)   =   (8.0)(l.0)   +    (10.0)(l0.0) 

=   8.0  +  100.0  =  108.0  Btu 

At  this  level  of  analysis  the  error  term  is  quite  large  and  dominated  by 
the  component  Ae'.  This  error  can  be  reduced  by  identifying  individual 
inputs  to  the  widget. 

Step  11(a).   The  analyst  obtains  additional  information  to 
classify  inputs  to  the  widget  by  sector.   Requirements  from  sectors  1 
and  3  are  found  to  be  typical  but  those  from  sector  2  are  felt  to  be 
atypical.   There  are  still  some  requirements  that  remain  unclassified. 
The  sector  input  values  and  uncertainties  are  found  to  be: 


Aa0  =  .1   Ae'  =  2.0  (Note:   Ae'  =  1.0) 


a1  =  2.0 

Aa     =   .1 

a2  =   3.9 

Aa2  =    .1 

a3  =   k.O 

Aa_  =    .8 

a     =  0.1 
o 

Aa     =    .1 
o 

New  estimates  of  energy  costs  and  first  order  uncertainties 
are  obtained  (inputs  ordered  in  counterclockwise  progression  from 
Fig.  h). 

Energy  Costs  =  lea     =   (U.0)(2.0)+(15.0)(U.0)+ 

(T.0)(3.9)  +  (8.0)(.l)  =  8.0  +  60.0  +  27.3  +  .8 
=  96.I  Btu 

Positive  Error  Term  =  E(eAa  +  Aea) 

=  [(U.o)(.i)  +  (.5)(2.o)]  +  [(15.0K.8)  +  (1.0MU.0)]  + 

[(7.0M.1)  +  (2.0)(3.9)]  +  [(8.0)(.l)  +  (IO.OH.1)] 

=  (.h   +  1.0)  +  (12.0  +  k.0)   +  (.7  +  7.8)  +  (.8  +  1.0) 
(sector  l)    (sector  3)    (sector  2)    (misc.) 

=  27.7  Btu 

The  error  term  is  nearly  a  fourth  of  the  one  obtained  in  the  first  step, 
but  is  still  about  25$  of  the  calculated  energy  cost.  The  major  component 
of  error  is  due  to  the  Aa  term  for  sector  3  inputs.  This  error  can  be 
reduced  by  returning  to  the  data  source  for  a  better  estimate  of  sector  3 
inputs. 

Step  11(b).   The  analyst  obtains  an  improved  sector  3  input 
specification.   Step  11(a)  is  repeated  except  that  Aa_  has  been  reduced 
from  .8  to  .1  (the  a_  estimate  is  found  to  be  the  same  as  before).   The 
resulting  energy  cost  and  uncertainty  are  calculated: 
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Energy  Costs  =  96. 1  Btu  (same  as  Step  11(a)) 

Positive  Error  Term  =  27.7  -  12.0  +  1.5  =  17.2  Btu 

old    new 


Now  that  the  12.0  error  term  is  reduced  to  1.5  the  major  contributor  to  the 
uncertainty  is  due  to  ktt    (the  7-8  term).   This  error  component  can  be 
reduced  by  specifying  the  inputs  required  to  produce  the  atypical 
sector  2  output. 


Step  III.   The  analyst  disaggregates  the  inputs  into  the 
atypical  sector  2  output.   The  following  typical  inputs  are  found: 

a   =  .1  Aa   =   .01 

ou  =  3.8  Aotp  =   .1 

Note  from  Fig.  5  that  there  are  now  five  inputs  crossing  the  system 
boundary.   Using  the  counterclockwise  accounting  scheme,  the  first  two 
and  last  one  are  identical  to  the  previous  step.   Incorporating  the 
new  information, 

Energy  Costs  =  lea     =  (U.0)(2.0)  +  (15.0)(U.0)  + 
(U.OK.1)  +  (7.0)(3.8)  +  (8.0K.1)  =  8.0  + 
60.0  +  .U  +  26.6  +  .8 
=  95.8  Btu 


-  19  - 

Positive  Error  Term     =     E(eAa  +  Aea) 

=  [(U.0)(.l)  +  (,5)(2.0)]  +  [(15.0)(.1)  +  (l.O)(U.O)]  + 
[(U.0)(.0l)  +  (.5)(.l)]  +  [(T.0)(.l)  +  (1.0)(3.8)]  + 
[(8.0)(.l)  +  (10.0)(.1)] 

=  (o.U  +  1.0)  +  (1.5  +  U.o)  +  (.oU  +  .05)  + 

(sector  l)    (sector  3)    (sector  l) 

(.7  +  3.8)  +  (.8  +  1.0) 
(sector  2)    (misc. ) 

=  13.29  Btu 

The  uncertainty  is  now  reduced  to  an  eighth  of  that  in  the  initial  step. 
The  major  components  of  error  are  due  to  Ae'  error  terms  for  the  sector  3 
inputs  (the  k.O   term)  and  the  (typical)  sector  2  input  into  the  atypical 
sector  2  requirement  (the  3.8  term).   Further  reduction  in  these  uncer- 
tainties can  only  result  from  an  improved  1-0  specification.  Therefore, 
this  seems  an  appropriate  point  to  terminate  the  analysis.   The  resulting 
energy  cost  of  the  widget  becomes: 

Energy  Cost   =  95.8  +  13.29  Btu. 

(in  an  actual  analysis  negative  uncertainties  would  be  treated  in  a 
manner  identical  to  the  positive  errors  considered  in  this  example.) 

The  table  below  summarizes  the  sample  case  described  above. 
Although  this  example  is  only  hypothetical,  note  that  the  reduction  in 
uncertainty  decreases  rapidly  as  the  analysis  progresses. 


-RO- 


TABLE 1.   Summary  of  three-sector  example. 


STEP       ACTION  TAKEN       ENERGY  COST    UNCERTAINTY   UNCERTAINTY 

(Btu)  (Btu)      REDUCTION  (Btu) 


80.0  108. 


Il(A)   Process  analysis  of 


widget  inputs 

96.1 

27.7 

81.3 

11(B) 

Tighten  Aa_ 

96.1 

17.2 

10.5 

III 

Process  analysis  of 
atypical  sector  2 

inputs 

95.8 

13.3 

3.9 

